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Abstract: Recent renormalization group calculations of the sparticle mass spectrum in 
the Minimal Supersymmetric Standard Model (MSSM) show that t — b—r Yukawa coupling 
unification at Mqut is possible when the mass spectra follow the pattern of a radiatively 
induced inverted scalar mass hierarchy. The calculation is entirely consistent with expec- 
tations from 50(10) SUSY GUT theories, with one exception: it seems to require MSSM 
Higgs soft term mass splitting at MquTj dubbed "just-so Higgs splitting" (HS) in the lit- 
erature, which apparently violates the 50(10) gauge symmetry. Here, we investigate three 
alternative effects: i). 50(10) D-term splitting, ii). inclusion of right hand neutrino in 
the RG calculation, and in), first/third generation scalar mass splitting. By combining 
all three effects (the DR3 model), we find t — b — t Yukawa unification at Mgut can be 
achieved at the 2.5% level. In the DR3 case, we expect lighter (and possibly detectable) 
third generation and heavy Higgs scalars than in the model with HS. In addition, the light 
bottom squark in DR3 should be dominantly a right state, while in the HS model, it is 
dominantly a left state. 
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1. Introduction 



The four LEP experiments performed precision measurements of the SU{3)c, SU{2)l and 
U{1)y gauge couphng constants of the Standard Model (SM) at energy scale Q = Mz [1]. 
It is an astonishing fact that the values of these couplings, evolved in energy from the weak 
scale to an energy scale Mqut — 2 x 10^^ GeV, nearly meet at a point under Minimal 
Supersymmetric Standard Model (MSSM) renormalization group (RG) evolution [2], while 
their unification fails badly under Standard Model RG evolution. This latter fact is often 
touted as indirect evidence that the MSSM (with weak-scale sparticle masses) is the correct 
effective field theory describing nature at energy scales between Mqut and M„eak) and 
further that nature may well be described by a supersymmetric grand unified theory (SUSY 
GUT) at energy scales above A/gut- 

While the gauge group SU{5) early on emerged as a leading GUT group candidate [3], 
the gauge group 50(10) appears to be much more compelling [4]. In SUSY S'O(IO) GUTs, 
a number of attractive features emerge. 

• All matter fields of a single generation are unified into the 16-dimensional spinor of 
50(10): thus, 50(10) unifies matter as well as forces. 

• The seemingly ad-hoc cancellation of triangle anomalies in the SM and in SU{5) is 
a simple mathematical fact in 50(10). 

• The 16th element of the 50(10) matter spinor is naturally occupied by a superfield 
iV"^ which contains a SM-gauge singlet right-handed neutrino (RHN) field. Upon 
breaking of 50(10), the RHN acquires a Major ana mass Mn which leads to the 
famous see-saw relation for neutrino masses [5]: nii, = {fvVuY I^N- 

• The fact that matter superfields lie in a spinor representation automatically leads 
to i?-parity conservation, since only superpotential terms of the form matter-matter- 
Higgs are allowed by the 50(10) symmetry, while the i?-parity violating matter- 
matter-matter or mxitter-Higgs products are not allowed. While 50(10) breaking 
may re-introduce i?-parity violation, many simple breaking schemes exactly preserve 
the i2-parity conserving structure. 

• 50(10) SUSY GUTs naturally explain why two Higgs doublets occur in nature at 
the weak scale. The 2 and 2* MSSM Higgs doublets lie in a 5 and 5* of 5C/(5), and 
the 10 of 50(10) naturally contains a 5 and 5* under restriction to SU{5). 

Thus, 50(10) SUSY GUTs provide an extremely compelling picture of what physics might 
look like around the GUT scale. 

Along with gauge coupling and matter unification, the simplest 50(10) SUSY GUT 
models also predict third generation t — h — t Yukawa coupling unification at the GUT 
scale [6]. To check this assertion, one must begin with the measured third generation 
fermion masses — mj, mi, and rrir — and calculate the associated Yukawa couplings at the 
weak scale. Then one may evolve the t — h — t Yukawa couplings up in energy to check 
whether they unify at MquTj just as the gauge couplings do. The values of the weak 
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scale Yukawa couplings depend strongly on the ratio of Higgs field vevs: tan/3 = Vu/vd- 
Furthermore, unlike the gauge couplings, the Yukawa couplings have a large dependence on 
sparticle loop corrections when passing from the SM to the MSSM effective field theories [7] . 
Thus, the t — b — r Yukawa coupling unification depends on the precise form of the sparticle 
mass spectra of the MSSM. This latter fact offers a consistency check: if sparticle masses 
are found with the expected pattern, it would be strongly suggestive that an 50(10) SUSY 
GUT model is valid around the GUT scale. 

Detailed calculations of when t — b — t Yukawa couplings unify at Mgut within the 
MSSM context have been performed by several groups [8-14]. They use an S'O(10)-inspired 
model parameter space given by 

mi6, mio, M|), mi/2, Aq, tan/3, sign(/x) (1.1) 

where mie is the GUT scale mass of all matter scalars, mio is the GUT scale mass of 
Higgs scalars, parametrizes Higgs mass splitting (HS) or possible scalar mass D-term 
splitting (DT) (the latter can arise from the breaking of S'O(IO) gauge symmetry), mx/2 
is the unified GUT scale gaugino mass, is the unified GUT scale soft SUSY breaking 
(SSB) trilinear term, tan /3 = Vu/vd is the weak scale ratio of Higgs field vevs, and is the 
superpotential Higgs bilinear term, whose magnitude — but not sign — is determined by the 
electroweak symmetry breaking (EWSB) scalar potential minimization conditions. 

In practice, the two Higgs field soft breaking terms — ^'h^ ^'h^ — cannot be 
degenerate at Mgut and still allow for an appropriate radiative breakdown of electroweak 
symmetry (REWSB) [15]. Effectively, rnjj^ must be less than rnjj^ at Mqut in order to 
give rnjj^ a head start in running towards negative values at Mweak- We parametrize the 
Higgs splitting as m|^^ ^ = mfg T 2M^ in accord with nomenclature for D-term splitting 
to scalar masses when a gauge symmetry undergoes a breaking which reduces the rank of 
the gauge group [16]. While the D-term splitting shoTild apply to matter scalar SSB terms 
as well, in practice, better Yukawa unification for /v, > is found when the splitting is 
restricted only to the Higgs SSB terms. The mass splitting applied only to Higgs scalars, 
and not to matter scalars, has been dubbed "just-so" Higgs splitting in the literature [10]. 

In previous work, the above parameter space was scanned over via random scans [8,11] 
and also by more efficient Markov Chain Monte Carlo (MCMC) scans [12] to search for 
Yukawa unified solutions using the Isasugra subprogram of Is A JET [17] for sparticle mass 
computations.^ The quantity 

j^^ maxUuhJr) (evaluated at Q = Mgut), (1-2) 
mm{ft, Jb, U) 

was examined, where solutions with i? ~ 1 gave apparent Yukawa coupling unification. 

For superpotential Higgs mass parameter > (as favored by (g — 2)^ measurements), 
Yukawa unified solutions with i? ~ 1 were found but only for special choices of GUT scale 
boundary conditions [8-12, 14, 19]: 

^Ref. [14] confirms the general structure of Yukawa- unified models also using the SoftSUSY [18] spectrum 
generator. 
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• mi6 3 - 15 TeV, 



• ~ -2mi6, 

• mio ~ 1.2mi6, 

• mi/2 < "116, 

• tan (3 ~ 50. 

Models with this sort of boundary conditions were derived even earher in the context of 
radiatively driven inverted scalar mass hierarchy models (RIMH) which attempt to recon- 
cile suppression of flavor-changing and CP-violating processes via a decoupling solution 
with naturalness via multi-TeV first/second generation and sub-TeV scale third generation 
scalars [20,21]. The Yukawa-unified spectral solutions were thus found in Refs. [11,12] to 
occur with the above peculiar choice of boundary conditions as long as mie was in the 
multi-ToV regime. 

Based on the above work [11,12], the sparticle mass spectra from Yukawa- unified SUSY 
models are characterized qualitatively by the following conditions: 

• first and second generation scalars have masses in the ~ 10 TeV regime, 

• third generation scalars, /i and ruA have masses in the few TeV regime (owing to the 
inverted scalar mass hierarchy), 

• the gluino has mass rrig ~ 300 — 500 GeV, 

• the lightest neutralino x° is nearly pure bino with mass typically m^o ~ 50 — 80 GeV. 

The presence of a bino-like Xi along with multi-TeV scalars gives rise to a neutralino 
cold dark matter (CDM) relic abundance that is typically in the range ^^oh'^ ~ 10^ — 
10^, i.e. far above [22] the WMAP measured [23] value VtcDuh'^ = 0.110 ± 0.006 by 
several orders of magnitude. A very compelling solution to the Yukawa-unified dark matter 
abundance problem occurs if one invokes the Peccei-Quinn solution [24] to the strong CP 
problem, which leads to dark matter being composed of an axion [25, 26]/axino [27,28] 
admixture [12,14,29], instead of neutralinos. The axino then can serve as the LSP instead 
of the lightest neutralino [30,31]. Cosmological solutions with a re-heat temperature Tr 
high enough to sustain non-thermal leptogenesis {Tr ~ 10^ — 10^ GeV) could most easily 
be found if the dominant component of the cold dark matter consisted of axions rather 
than axinos. 

The above calculational results show that Yukawa unified solutions are compatible with 
the MSSM as the effective theory between the GUT and weaJi mass scales, but for a very 
constrained form of the sparticle mass spectrum. Indeed, the entire scheme is compatible 
with expectations from an 50(10) SUSY GUT theory with GUT symmetry broken at 
the scale MquTj save for one feature. Naively, 5*0(10) symmetry implies the two Higgs 
masses should be degenerate at Mqut, at tree-level: i.e. rn?fj^ = "^l^^- The mechanism 
for the large GUT scale Higgs soft term splitting is unknown, and violates the 50(10) 
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symmetry. There is of course a well-known source of Higgs soft term splitting: namely, 
the D-term contribution to scalar masses which is induced when the rank of the gauge 
group is reduced from 5 of 50(10) to 4 of SU{5) or the SM. The D-term contributions are 
necessarily proportional to the charge of the subgroup U{l)x, and are to be included in 
all scalar fields carrying a U{l)x charge. At Q = Mqut, the gauge symmmetry breaking 
induces a scalar mass contribution 



where the contribution M|, is effectively a free parameter of order the weak scale, and 
whose value depends on the details of 5*0(10) breaking. It can take positive or negative 
values.^ It was found in Refs. [9-11] that calculationally, the spectral solutions with the 
best t — b — T Yukawa unification occurred when the L>-term splitting was applied only to 
the Higgs scalar s, and not to the other matter scalars. 

In this paper, we re-visit the question of the HS case versus the DT splitting case in 
t — b — T Yukawa-unified models. We search for Yukawa-unified solutions while including 
several effects which are all consistent with the general framework of simple 50(10) SUSY 
GUT models: i.) application of full DT splitting to all scalar masses, ii.) inclusion of 
neutrino Yukawa coupling effects (RHN) [10, 21], and Hi.) inclusion of mass splitting 
between the third generation, versus the first two generations of matter scalars (3GS). We 
scan over 50(10) model parameter space using the MCMC technique, which provides an 
efficient search for the best Yukawa unified solutions. We find that each of the above three 
effects acts to improve the degree of Yukawa unification compared to results without the 
effects, but none of them work as well as the just-so HS model. However, using all three 
effects simultaneously (the DR3 model) does allow us to reach Yukawa-unified solutions 
with R ~ 1.025, i.e. Yukawa unification to below the 3% level. The remaining last few 
percent might then be compatible with expected GUT scale threshold effects (which are 
of course model dependent) and intrinsic theory error in our 2-loop RGE calculations. We 
find that the sTipcrparticlc mass spectra using the DR3 model is somewhat modified from 
solutions using just-so HS. In particular, for a given value of mig, the predicted value of 

and m^^ are much lighter than in the just-so HS prediction. For mig ~ 10 TeV, niA 
and can extend down to or even below the 1 TeV level, and may be accessible to 
LHC searches and/or to searches at a future CERN e+e~ Linear Collider (CLIC), where 
center-of-mass energies of order -s/s ~ 3 — 5 TeV are proposed. The spectral differences, 
and the expected left-right composition of the lightest bottom squark, may offer a means 
to distinguish between the just-so HS model and the DR3 model. 



^Sum rules for sparticle masses as a test of the underlying 50(10) are discussed in [32]. 



rn-D = 'mi 




mle + Ml 
mle - 3M|, 



(1.3) 
(1.4) 
(1.5) 
(1.6) 
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2. Higgs mass splitting and radiative EWSB in Yukawa unified models 



2.1 Radiative EWSB in Yukawa unified models 

MSSM models with t — b — t Yukawa coupling unification and degenerate (Higgs) soft 
masses at Q = Mqut face a well-known problem when one attempts to generate realistic 
sparticle mass spectra: electroweak symmetry fails to be appropriately broken [6,15]. The 
problem can be seen by examining the one-loop Higgs soft mass RGEs: 

-^glMf - SglMi - ^gjS + Zf^X, + /^X,) , (2.1) 

-^glMl - SglMi + ^^^2^ + 3f^Xt] , (2.2) 





2 


dt 


167r2 


dm?TT 


2 


dt 


167r2 



where 



Xt=ml,^ + m}^+ml^+A], (2.3) 
Xb = + m?^ + ml^ + Al , (2.4) 
Xr = ml^ + m?^ + ml^ + Al , (2.5) 



and 



i^Hn - "i/^d + [^Q -^L- 2m^ + + m|] . (2.6) 

The right-hand side terms with negative co-efficients give an upwards push to the Higgs soft 
masses during evolution from Mqut to M^eak) while the positive terms give a downwards 
push. 

At the weak scale, where the Higgs effective potential is minimized, the EWSB mini- 
mization conditions require that (at tree-level) 

^ (m^^+mf,^ + 2/x2)sin2/3 ^ ^^^^^ 

" - (tan2/3-l) ^^■^> 

The first of these determines the weak scale value of B in terms of tan the second relation 
determines the magnitude, but not the sign, of the superpotential /x term. At moderate- 
to-large tan /? values and ^ M'^, the second relation also gives approximately fj,'^ ~ 
, and we see that rn^j^ must be driven to negative values to accommodate successful 
REWSB. In models with Yukawa couplings ft > fb, fr, the m^j^ term is pushed to negative 
values by the large value of ft in Eq. ( |2.2D , resulting in successful EWSB. In contrast, in 
models with t — b — r Yukawa unification, the Yukawa coupling terms on the right-hand side 
of the fn^j^ equation are larger than the corresponding terms in the rn?ij^ equation, resulting 
in rn^H^ being driven more negative than tn^u^ at the weak scale. If w-^^ < rn?^^ i&v? (3, then 
IJ? < 0, signaling an inappropriate EWSB. The solution to this dilemma so far in Yukawa 
unified models is to provide the rn'jj^ term a head-start in running to negative values at 
the weak scale by adopting Higgs splitting such that rnj^^ < at the GUT scale. 
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While just-so HS applies a splitting (that violates 50(10) gauge symmetry) only to 
Higgs soft masses, and leaves the remaining GUT scale scalar masses fixed at m-ig, the 
expected splitting due to D-terms arising from 50(10) breaking at Mqut apply to matter 
scalars as well as Higgs scalars as given in Eqs. (|L3|)-([L^). The problem arising from DT 
splitting is that the rnj^ terms are substantially reduced already at MguT) and can get 
driven tachyonic at M^^ak through the RIMH mechanism. What is needed in terms of DT 
splitting is a large enough Higgs splitting to facilitate REWSB, but not so large a splitting 
that m? is driven tachyonic at M^eak- 

OR 

2.2 The role of neutrino Yukawa couplings 

An alternative mechanism to aide iTi'jj^ being driven negative is to balance the right-hand- 
side Yukawa push in Eqs. ( |2.lD -(^^) by incorporating the effect of the third generation 
neutrino Yukawa coupling. In this case, the "i^^ soft term RGE is modified to 



drn^T 



dt 167r2 



5 



(2. 



with Xi, = mj^ + fnf,^ + fnjj^ + A^, and the RGE for is unchanged. We see that 
for unified Yukawa couplings, now Tn'j^^ receives an additional downward push from the f,^ 
term [10,21]. The fu term contributes to the running at energy scales Mn < Q < Mgut- 
For Q < Mj\f, the RHN states are integrated out of the effective theory, and the RGEs 
revert to the MSSM form of Eq. ( ^ . 

2.3 Third generation splitting 

A third effect that can improve the implementation of EWSB in Yukawa unified models 
is to allow for non-degeneracy in first /second versus third generation scalar masses: third 
generation mass splitting (3GS). In general, generational non-degeneracy can lead to flavor 
violation in excess of experimental bounds [33]. The bounds on flavor-changing neutral 
currents (FCNC) apply most strongly to splitting amongst the flrst and second generation 
soft terms; contraints on third generation FCNCs are much less restrictive [34]. Here, 
we will maintain degeneracy between mi6(l) and mi6(2), although some small breaking 
in these terms is allowed, especially in our case where scalar masses are quite heavy and 
typically in the 10 TeV range, where we also have FCNC suppression via decoupling. Here, 
we will adopt rriiQ^l) = mi6(2) ^ mi6(3). 

At tree level, we would expect 3GS to not affect Yukawa coupling evolution and EWSB, 
since the Higgs sector only couples strongly to third generation scalars, and is independent 
of flrst/second generation scalar masses. However, at two-loop level, the scalar mass RGEs 
have the form given by [35] 

dt ~ 167r2^-?+ (167r2)2^-?' ^^'^^^ 

where t = InQ, i = Qj, Uj, Dj, Lj and Ej, and j = 1 — 3 is a generation index. Two 
loop terms are formally suppressed relative to one loop terms by the square of a coupling 
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constant as well as an additional loop factor of IGvr^. However, these two loop terms include 
contributions from all scalars. Specifically, the two loop /3 functions include, 

P3 ^ "'i93^3 + ^92(^2 + Cigfai, (2.11) 

where 

(71 = ^gl{3iml^ + m^J + Tr[uil + 3mi + Sm^, + 2inl + 6m|]}, 

c^2 = gl{m\j^ + mjj^ + Tr[3m^Q + mi]}, and 
0-3 = fiffTr [2m^ + niii + m|,], 

and the m? are squared mass matrices in generation space. The numerical coefficients 
Oj, bi and Cj are related to the quantum numbers of the scalar fields, but are all positive 
quantities. Incorporation of large, multi-TeV masses for the first and second generation 
scalars leads to an overall positive, non-negligible contribution to the slope of SSB mass 
trajectories versus energy scale [36]. Although formally a two loop effect, the smallness 
of the couplings is compensated by the multi-TeV scale values of masses for the first two 
generations of scalars. In running from Mqut to M^eak; the two- loop terms result in an 
overall reduction of scalar masses, and its effect depends on the quantum numbers of the 
various scalar fields. 

Generational non-degeneracy of scalar masses, especially for the 3GS scenario, is nat- 
ural in SO{10) SUSY GUT models where above-the-GUT-scale-running is allowed. In this 
case, the unified third generation Yukawa coupling acts to suppress mi6(3) with respect 
to m-i6(l,2), even if the three generations are degenerate at some higher scale, e.g. at the 
Planck scale Mp. 

SO{W) model RGEs are presented in Sec. 6.1 of Ref. [21]. As an example, we show in 
Fig. |l]the evolution of mfgi ^-16(1) 2) and mi6(3) from Mp to Mqut for model parameters 
as depicted in the caption. We see that a splitting of order 25% is possible at Mqut for 
parameter choices as befitting Yukawa unified models. The natural splitting here is then 
that mi6(3) < mi6(l,2), assuming degeneracy at Mp. In our following calculations, we 
will merely implement mi6(3) as a free parameter as distinct from mi6(l,2). 

3. Numerical results for Yukawa unified models 
3.1 MCMC scan of parameter space 

In this section, in addition to the just-so HS model, we will also scan over the SO{W) 
model parameter space as given by 

mi6(l,2), mi6(3), mio, Md, mi/2, ^o, tan/?, sign{n), (3.1) 

where we also update our boundary conditions of mb{Mz)^^ = 2.416, as{Mz) = 0.118 
and mt = 173.1 GeV; the latter in accord with recent Tevatron measurements [37]. In 
addition, we augment the parameter space with the RHN parameter set 

Mn^s, fur, ^ur, m^m (3-2) 
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Figure 1: Renormalization group evolution of soft SUSY breaking terms mi6(l,2), mi6(3) and 
TOio from Mp to Mqut in a simple S'O(IO) SUSY GUT. We take gcuT = 0.682, four = 0.56, 
TO16 = 10 TeV, mio = 12.05 TeV and Aq = -19.947 TeV (point B of Table 2 of Ref. [14]). 



where we assume fu = ft^^Q = -^guTj = and ruj^^^ is as given in Eq. (L5) by DT 
splitting. For M^fg , we are guided by the simple see-saw relation 



where for tan/5 ~ 50, we have Wu(msusY) 



171.6 GeV, 



Vdimsvsy) 



(3.3) 



3.5 GeV and 



/^^ ~ 0.54. Then, cosmological bounds on the sum of neutrino masses implies ^ ~ 1 



eV. This imphes Mjv, ~ 10^^ GeV. We will adopt here Mjvg = 10^^ GeV to maximize 
the downward push of f^^ consistent with bounds on neutrino masses. The superparticle 
mass spectrum is then generated using Isajet 7.79 [17], which includes full two-loop RGE 
running, implementation of the RG-improved 1-loop effective potential for EWSB, and 
full 1-loop corrections to all sparticle masses. We scan over the parameter space using 
the MCMC method with a Metropolis sampling algorithm as in [12], which provides an 
optimized search for parameter space points with the lowest R values. 

In the scans, we require that the mass limits from direct SUSY [38] and Higgs [39] 
searches at LEP be observed (additional limits from Tevatron searches do not affect the 
solutions with small R). Moreover, we take into account the constraints from the branching 
fractions for the 6 — > 57 and Bg decays. The measured branching ratio of the 

inclusive radiative B decay is BR(6 57) = (3.52 ± 0.23 ± 0.09) x lO"'' [40], and the SM 
theoretical prediction BR(6 — > 57)^^ = (3.15 it 0.23) x 10~^ [41]. Combining experimental 
and theoretical errors in quadrature, we take 2.85 < BR(6 — > 57) x 10^ < 4.19 at 2a together 
with the 95% CL upper limit BR(Ss ^ < 5.8 x 10"^ from CDF [42]. We adopt the 

Isajet Isatools [43] program for the calculations of BR(6 sj) and BR(i?s fi^n^). 

Our results for parameter space scans with minimized R values are summarized in Table 
|l[ We first see that t — h — r Yukawa unification is not possible in the mSUGRA model: in 
this case, Rmin = 1-35. However, adopting 50(10) parameter space with RHN parameters 
(but with Md = 0) actually allows considerable improvement with R values as low as 1.19 to 
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model Rmin 

mSUGRA 1.35 

mSUGRA+RHN 1.19 

5O(10)+RHN, Md = 1.19 

5O(10)+DT 1.08 

5O(10)+3GS 1.30 

5O(10)+DT+3GS 1.06 

5O(10)+DT+RHN 1.04 

5O(10)+RHN+3GS 1.17 

5O(10)+DT+RHN+3GS (DR3) 1.025 

5O(10)+HS 1.0 



Table 1: Minimal R values obtained from MCMC scans for < — 6 — t Yukawa unification in various 
model parameter space choices. 

be attained. The same Rmin can be obtained in mSTJGRA-t-RHN. Allowing 5'O(10)-t-DT 
splitting (but no RHN), we find Rmin = 1-08, while S'O(10)+HS yields Rmin = 1-0, as 
noted in many previous studies. An allowance of DT+3GS yields a value Rmin = 1-06 
while allowing DT+RHN gives Rmin = 1-04. S'O(IO) scans with DT+RHN+3GS (the 
DR3 model) allow us to generate models with Rmin values down to 1.025. The remaining 
few per cent may be accounted for either by GUT scale threshold corrections, or theoretical 
error due to our imperfect modeling of Yukawa coupling boundary conditions and evolution 
over 13 — 14 orders of magnitude. If these combined effects are at the level of few per cent, 
then the 50(10) model with D-term splitting plus RHN and/or 3GS can be seen to be 
in accord with Yukawa unified models for > 0. In the following, we concentrate on 
the DR3 model, which combines DT+RHN+3GS and gives the best prospects for Yukawa 
unification, and compare it to the "just-so" Higgs splitting. 

Next, we examine which parameter choices lead to t — 6— r Yukawa coupling unification 
in the DR3 model. To this end, we plot the locus of Yukawa unified solutions with R < 1.05 
in Fig. H in the mi6(3)/mi6 vs. Md plane. ^ Here and in the following, mig = m-i6(l,2) 
for simplicity. The red points indicate results from the just-so HS model. They necessarily 
all have mi6(3)/mi6 = 1, and so form a vertical line through the plot. The values of 
in the HS model actually range up to ~ 5 TeV. In contrast, we find that the DR3 model 
yields Yukawa-unified solutions provided that miQ{3) ~ (0.8 — 1.05)mi6(l, 2): i.e. the 
third generation scalars are favored to be at somewhat lower GUT scale masses than their 
first /second generation counterparts. In addition, we see that Md is now restricted to lie 
in the 1 — 2 TeV range: significantly less than the splitting needed for the HS model. 

In Fig. Ill, we show the locus of points with R < 1.05 in the Md vs. uiiq plane. Here 
we see that for the HS model, values of Md ~ 0.33mi6 are needed. In the case of the DR3 
model, we see that the value of Md needed also grows with mig. But in this case, we find 
instead that Md ~ 0.13mi6. The amount of Higgs splitting needed in the DR3 model is 
much less than in the HS model since additional Higgs splitting comes from the effect of 

^Note that in our case is always positive since we need m^i^ < at the GUT scale. 
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Figure 2: Yukawa unified solutions with R < 1.05 from "just-so" HS model (red points) and the 
DR3 model (blue points) in the mi6(3)/TOi6 vs. Md plane. Light blue/light red points are excluded 
by B-physics constraints. Note that Md extends up to about 5 TeV in the HS case, cf. Fig. ||. 



the RHN and 3GS. Other than the relative reduction in Md needed for a given value of 
771-16(1,2), the usual SSB mass relations for Yukawa-unified HS model still hold. 

Note also that in the DR3 model, Yukawa-unified solutions occur only for mi6(l, 2) ~ 
8 TeV, while in the HS model they also occur for smaller values of mig- However, the HS 
points with mig ~ 4 — 8 TeV are almost all excluded by BR(6 57) and/or BR(i?s — > 
fi^fi^). This observation is also made by Altmannshofer et al. [13]. The DR3 model points, 
on the other hand, are much less affected by the i3-physics constraints. We will come back 
to this in section 3.3. 

3.2 SUSY particle mass spectrum 

Given that both the HS and DR3 models lead to t — b — t unification at MquTj the 
next question is whether it is possible to physically distinguish between these models at 
experiments. Several differences between the SUSY particle mass spectrum lead us to 
believe that the models are at least in principle distinguishable. The first point is that — in 
the case of the DR3 model — the GUT scale soft masses for rnjj and m| are diminished by 
— 3M|) relative to the value of miQ. The first and second generation values oirriD and niL 
are expected to be in the multi-TeV regime, and so their mass diminution by D-terms isn't 
likely to be visible at any collider expected to operate in the near future. However, the 
third generation scalar masses are driven to weak scale values by the RIMH mechanism, 
and are expected to be in the 1-2 TeV regime. Thus, we would expect the third generation 
hn and fi masses to be diminished with respect to expectations from the HS model. This 
effect should be most noticable in the bottom squark sector, since in the tau slepton sector. 
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Figure 3: Yukawa unified solutions with R < 1.05 from "just-so" HS model (red points) and 
models with DR3 splitting (blue points) in the Md vs. mig plane; points excluded by _B-physics 
constraints are shown in lighter colour. 



we usually expect (based on the form of the MSSM RGEs) m^^ > ?Tif^, whereas in the 
sbottom sector, we expect m^^ < m^^ . 

In Fig. ^, we show the value of m^^ vs. rriiQ for Yukawa-unified models with R < 1.05 
in the HS (red dots) and DR3 (blue dots) cases. Points which obey the mass limits but 
do not pass i?-physics constraints are again shown in lighter colour. We see that for a 
given value of mie, the value of m^^ is smaller in the DR3 case than in the HS case. 
Naively, one might expect the opposite result, since in the DR3 model both m? and m? 
are increased by However, two effects act counter to the D-term. First, there is the 

third generation mass splitting, which typically reduces mi6(3) relative to miQ{l). Second, 
there is an additional contribution to RG running of top squark soft masses from the S- 
term, Eq. ( p.6D , which is zero for models with strict universality, but which is non-zero 
for models with Higgs mass splitting. For our case with ni'jj^ < fn^j^, the S-term gives 
an upwards push to the top squark soft mass evolution. Since the Higgs splitting is much 
less in the DR3 model, there is a smaller upwards push from the S'-term, and this effect 
coupled with 3GS wins out over the increased mass due to the D-term, thus giving the 
DR3 models typically a smaller ti mass than in the HS case. We note here that while the 
value of is smaller in the DR3 case compared to HS- for a given value of mig- the 
value of mi6 will not be easily measureable, and so the mj^/mig ratio is not likely a good 
discriminator between the two models: for instance, if a value of m^^ ~ 1.5 TeV is found 
at some future experiment, it will be difficult to determine if it is consistent with the HS 
model with mie ~ 7 TeV, or with the DR3 model with mig ~ 9 TeV. 



- 11 - 



4.5 

4 

3.5 

_ 3 
> 

^2.5 

1.5; 
1; 

0.5r 



DR3, allowed by B obs 
HS, allowed by B obs 
DR3, excluded by B obs 
HS, excluded by B obs 







t ,1 1 -'" 




I I I I I I I I 







% 4 5 6 7 8 9 10 11 12 13 
m^6(TeV) 

Figure 4: Yukawa unified solutions with R < 1.05 from "just-so" HS model (red points) and 
models with DR3 splitting (blue points) in the mig vs. rrif^plane; points excluded by _B-physics 
constraints are shown in lighter colour. 



In Fig. ^, we plot the value of m^^ expected from Yukawa-unified models with HS (red 
dots) and DR3 (blue dots), where again we only plot solutions with R < 1.05. We see that 
in the HS case, for a given value of mie, the value of m^^^ is always lowest in the DR3 case. 
In fact, for mig ~ 10 TeV, we expect m^^ ~ 1 — 2 TeV, while in the HS case, m^^ ~ 3 — 4 
TeV. 

The mixing angle of bi is also strongly affected. Here, following the notation of Ref. [44], 
we have bi = cosdJ}L~s,\ii9j)R. As an illustrative example, we list in Table |2| the spectrum 
from a HS Yukawa-unified point with mig = 10 TeV, and a DR3 model case with mig = 
11.8 TeV. The evolution of all four Yukawa coupings are shown for the DR3 point in Fig. |6|. 

In the HS case, 61 is ~ 10% while in the DR3 case, hi is 99.8% This comes 
from the fact that the D-term mass contribution pushes m^^ up, and down. We also 
show in Fig. the 6-squark mixing angle Ob versus fn^_^. In this plot, we see the value of 
6ij ~ 1.5 in the DR3 case, which means the 61 is dominantly Meanwhile, the red points 
indicate that in the HS model, the 61 is dominantly b^. 

The composition of 61 in principle might be measureable at the LHC if a bottom squark 
production event sample can be isolated, and the 61 branching fractions can be measured. 
Here note that the bi decays into wino-like charginos and neutralinos and/or into Wti, 
while the 6r does not. To give a concrete example, we compare the 61 decays at the DR3 
point of Table § to those of a HS point with a very similar 61 mass: point A of Ref. [14]. 
This latter point has mie = 5 TeV, m^^ = 1322 GeV, m^^ = 834 GeV, rrig = 363 GeV, 

= 109 GeV and m^o = 50 GeV. Several relevant 61 branching fractions are listed 

Xi 5X2 xi 
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Figure 5: Yukawa unified solutions with R < 1.05 from the "just-so" HS model (red points) and 
models with DR3 splitting (blue points) in the mig vs. plane; points excluded by _B-physics 
constraints are shown in lighter colour. 




in Table ^ As can be seen, in the DR3 model, 61 decays nearly 100% of the time into bg, 
while in the HS model, there is a sizable branching into txi and Wti states. 

While it is conceivable that the L — R composition of bi might be measured at LHC, 
the measurement would likely be very difficult and intricate. However, a measurement of 
the composition of 61 would likely be quite straightforward at a linear e~^e~ collider with 
sufficient energy to produce bibi pairs. First, the branching fractions of hi would likely be 
much easier to dis-entangle at an e~^e~ collider than at the LHC. Second, a linear e~^e~ 
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collider is likely to be constructed with polarizable electron beams. The total e~^e~ bibi 
cross section will be very sensitive to the beam polarization, and the composition of the 
bi. The beam polarization-dependent sparticle pair production cross sections have been 
calculated in Ref. [45], and the results are plotted versus the beam polarization Pi(e~) in 
Fig. |8[ Here, we see that a{e~^e~ — > bibi) peaks at PL{e~) = +1 in the case of the HS 
model, whereas just the opposite peak at Pi{e~) = —1 is expected in the DR3 case. 

In Fig. ^, we show the value of uia vs. uiiq for Yukawa-unified models with R < 1.05 
in the HS and DR3 cases. For moderate to large tan /3, at tree level we roughly expect 
that ~ '^Hd ~ ^Hu^ ^'^^ that the value of m\ is nearly equal to the weak scale mass 
splitting between the two Higgs soft masses. In the DR3 model, where a much smaller 
Higgs splitting is needed at the GUT scale in order to accomplish REWSB, we also find a 
significantly smaller value of niA expected for a given value of miQ. For miQ ~ 10 TeV, we 
typically get uia ~ 1 TeV in the DR3 model, while rriA ~ 3 TeV in the HS model. 

The value of niA may be readily established at the LHC, especially for the large 
tan/3 ~ 50 case expected for Yukawa-unified models. In the large tan/3 case, the 6-quark 
Yukawa coupling is large, and production cross sections for A are enhanced, both via glue- 
glue fusion (triangle diagrams) and via bA and bbA production. Then the A is typically 
expected to decay into modes such as bb, t^t^ and fJ-^fJ-^- The first two of these offer a 
rough mass bump with which to reconstruct rriA', the latter mode into ^'^ suffers from 
a small branching fraction, but may offer a sharper mass bump since all the decay prod- 
ucts are easily detected [46]. However, for extremely high energy muons, the momentum 
resolution — determined by track bending in the detector magnetic field — gets more dif- 
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parameter 


Pt. B [14J 


DR3 


mi6(l,2) 


10000 


11805.6 


mi6(3) 


10000 


10840.1 


mio 


12053.5 


13903.3 


Md 


3287.1 


1850.6 




43.9442 


27.414 


A 

^0 


-19947.3 


-22786.2 


tan (3 


50.398 


50.002 


R 


1.025 


1.027 




3132.6 


2183.4 


m~g 


351. z 


ozl.4 




QQ79 1 


1 1 Q1 zL 9 


"Hi 


2756.5 


2421.6 




3377.1 


1359.5 


HI ~ 


10094.7 


1 1 968 5 


xf 


116.4 


114.5 


X2 


113.8 


114.2 




49.2 


46.5 


ruA 


1825.9 


668.3 


ruh 


127.8 


128.6 


Oh (radians) 


0.329 


1.53 



Table 2: Masses in GeV units and parameters for Yukawa-unified point B of Ref. [14] with just-so 
HS, and a point with the DR3 model using Mjv = 10^^ GeV. We also give the 6-squark mixing 
angle. 



parameter Pt. A [14] DR3 







1321.8 


1359.5 


bi - 


' tXi 


8% 


0.1% 


h - 


> Wil 


30% 




h - 


>h~g 


55% 


99% 



Table 3: Sbottom mass and branching fractions for a HS model (point A of Ref. [14]) and for the 
DR3 point from Table 2. Note that the two points have very similar m^^, but the hi is mainly a 
left-squark in the HS case, while it is mainly a right-squark in the DR3 case. 



ficult at high energies. Evaluation of the A mass bump in all three of these modes may 
allow good resolution on the A mass reconstruction. If a measurement of ^ — > ^i'^ ^~ is 
possible, then it may also be possible to extract information on the width Ta, which is 
very sensitive to the value of tan (5. 



3.3 Predictions for i?-physics observables 



For completeness, we present in Fig. 10 the explicit results for the branching ratios of the 
6 — > S7 and Bg — > /i decays for the Yukawa-unified solutions with R < 1.05 of the 
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Figure 8: Sbottom production cross section, a{e^e~ — > feifci), in fb versus electron beam polar- 
ization Pi(e~) for the CLIC accelerator with ^/s = 3 TeV. The blue curve corresponds to a DR3 
model point while the red curve corresponds to the HS model point A of Ref. [14]. 
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Figure 9: Yukawa unified solutions with R < 1.05 from "just-so" HS model (red points) and 
models with DR3 splitting (blue points) in the mig vs. tua plane; points excluded by _B-physics 
constraints are shown in lighter colour. 



previous section. 

For the HS model, we see that most of the solutions with mie ~ 8 TeV lead to too 
low a value of BR(6 — > 57), and hence are excluded. Of the remaining HS points with 
mi6 ~ 8 TeV, a large fraction has too high a BR(i3s fi~^fi^). In the end, only a 
few points with mig ~ 8 TeV survive and they have BR(i?s — ^ A'^A^ ) = O(10~^). The 
HS points with miQ > 8 TeV mostly comply with the i?-physics constraints and have 
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Figure 10: Predictions for BR(6 57) vs. mie (left plot) and BK{Bs — > fi^fj.~) vs. rriiQ (right 
plot) from Yukawa-unified models in both the HS (red dots) and DR3 (blue dots) models. Also 
shown are the 2a limits used in our analysis. 



BR{Bs ^+Ai-) = 0(10-*^). 

The situation is different for the DR3 model points, for which mie is always greater 
than 8 TeV. These points cluster around BR(6 — > 57) ~ 3 x 10~^, with only a very small 
fraction (less than 1%) excluded by a too high BR(i?s — > fi^fi^). Indeed the excluded 
points are those with the lowest values, cf. Fig. |^. Interestingly, owing to the lower 
niA, the DR3 model predicts BR(i?s —>■ iJ-^fi~) = 0(10^^), which may well be probed in 
the near future at the Fermilab Tevatron collider. 



4. Summary and conclusions 

In this paper, we have re- investigated the issue of electroweak symmetry breaking in 
Yukawa-unified SUSY models. It is well known that perfect Yukawa coupling unifica- 
tion can be achieved in simple SO (10) SUSY GUTs with universal soft-breaking terms 
for gauginos and sfermions, but non-universal Higgs-mass terms. However, it is hard to 
understand how this "just-so" Higgs splitting may arise only in the Higgs multiplets, and 
not in the corrsponding matter multiplets, as is expected in D-term splitting. We have 
found that Yukawa coupling unification good to few per cent can be achieved in the case 
of D-term splitting but only if one also allows for the presence of the neutrino Yukawa 
coupling, along with first /second versus third generation matter scalar splitting. Each of 
these three features — L>-term splitting, right hand neutrino Yukawa couplings and third 
generation splitting, together comprising the DR3 model — are to be expected in simple 
GUT models based on 5*0(10). In this case, the DR3 model may be considered more 
satisfying from the 50(10) GUT point of view than the HS model. 

The two models lead to many similarities in the expected sparticle mass spectra, but 
also some important differences. Regarding similarities, both lead to a split spectra with 
first/second generation scalars in the ~ 10 TeV regime, with third generation and heavy 
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Higgs scalars in the few TeV regime, along with a very Ught spectrum of gauginos. In 
particular, with gluinos expected in the 300 — 500 GeV mass range, a robust variety of 
gluino pair production events are expected at the CERN LHC [47]. Regarding differences, 
the amount of Higgs splitting for a given value of mig is expected to be much less in the 
DR3 model, leading to much lighter values of itia, mn and mjj±. These heavier Higgs 
states stand a much higher chance to be detectable at LHC in the DR3 model, as compard 
to the HS model. In addition, the lightest bottom squark is expected to be much lighter for 
a given value of mig in the DR3 model, as compared to the HS model. Also, the smoking 
gun difference is that the 6i should be predominantly a right-squark in DR3, while it is 
expected to be dominantly a left-squark in the HS model. This can in principle be detected 
at LHC due to the different bi branching fractions which are expected; however, in practise, 
this differentiation is likely to be a difficult enterprise. It will be much simpler at a CLIC- 
type e~^e~ linear collider operating with ^/s > 2m^^. In this case, the bi production and 

decay modes should be more readily identified, and especially the total bibi production 
cross section will depend on electron beam polarization in very different fashions for the 
two models. In this case, the two models should then be easily distinguishable. 

We also evaluated predictions for BR(6 — > 57) and BR(i?s in the HS and DR3 

models. In the HS model, most points with rriiQ < 8 TeV are excluded by these constraints. 
In the case of the DR3 model, which apparently requires mie ~ 8 TeV, 99% of the points 
are allowed by the i?-decay constraints. But the predicted rate for the Bg IJ'~^IJ'~ decay 
is just below its current experimental limit, and may well be probed in the near future as 
data continues to accrue at the Fermilab Tevatron collider. 

Finally we note that a large RHN Yukawa coupling, as assumed in this study, can 
lead to lepton flavor violation (LFV) by generating off-diagonal LFV terms in the charged- 
slepton mass matrix through RG running [48,49]. This can lead to additional constraints 
which may further sharpen the predictions for Yukawa-unified models. This issue is left 
for future work [50]. 

Note added: As this paper was being finalized, some related papers appeared on 
Yukawa coupling unification in models without universality in gaugino masses [51] and 
^-terms [52]. 
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